function plot_loca(loc_BS_ATG,loc_airplane,loc_BS,loc_user,num_terres_BS_layer,dist_BS_half)
range_point = range_point_generate(num_terres_BS_layer,dist_BS_half);
%
LineWidth = 1.5;
markersize = 10;
% 4DBEEE
ff=figure;
plot(loc_BS_ATG(:,1),loc_BS_ATG(:,2),'p','color','#D95319','Markersize',markersize,'LineWidth',LineWidth);hold on
plot(loc_airplane(:,1),loc_airplane(:,2),'<--','color','#4DBEEE','Markersize',markersize,'LineWidth',LineWidth);hold on
% plot(loc_BS_all(:,1),loc_BS_all(:,2),'*','color','#0072BD','Markersize',markersize,'LineWidth',LineWidth);hold on
plot(loc_BS(:,1),loc_BS(:,2),'*','color','#0072BD','Markersize',markersize,'LineWidth',LineWidth);hold on
plot(loc_user(:,1),loc_user(:,2),'x','color','#7E2F8E','Markersize',5,'LineWidth',LineWidth);hold on
plot(range_point(:,1), range_point(:,2),'--','color','#77AC30','linewidth',LineWidth);hold on;
legend('ATG BS','airplane','Terrestrial BS','Terrestrial user',...
    'Considered area','FontSize',10,'Location','southeast');
xlabel('x (km)')
ylabel('y (km)')
% set(gca,'xtick',[min(xlabel_set_1):max(xlabel_set_1)]);
% set(gca,'ytick',0:20);
grid on
axis equal
% ylim([-8.5*dist_BS_half,8.5*dist_BS_half]);
ylim([-6*dist_BS_half,6*dist_BS_half]);
set(gca, 'FontName', 'Times New Roman','FontSize',12)
% set(ff,'Position',[1000,394,500,400]);
end

function range_point = range_point_generate(num_terres_BS_layer,dist_BS_half)
range_point = [pi/3,num_terres_BS_layer*2+1;...
    pi/3*2,num_terres_BS_layer*2+1;...
    pi/3*3,num_terres_BS_layer*2+1;...
    pi/3*4,num_terres_BS_layer*2+1;...
    pi/3*5,num_terres_BS_layer*2+1;...
    pi/3*6,num_terres_BS_layer*2+1;...
    pi/3,num_terres_BS_layer*2+1];
[xxx,yyy]=pol2cart(range_point(:,1)',range_point(:,2)'*dist_BS_half);
range_point(:,1) = xxx'+num_terres_BS_layer*2*dist_BS_half;
range_point(:,2) = yyy';
% index_plot = [1,2;2,3;3,4;4,5;5,6;6,1];
end